Fractional quantum Hall states of few bosonic atoms 
in geometric gauge fields 



B. Julia-Diaz 1,2 , T. Grafi 2 , N. Barberan 1 , and M. Lewenstein 2 ' 3 

1 Dcpartamcnt d'Estructura i Constituents de la Materia, Facultat de Fisica, 
Universitat de Barcelona, E-08028, Barcelona, Spain 

2 ICFO-Institut de Ciencies Fotoniques, Pare Mediterrani de la Tecnologia, 08860 
Barcelona, Spain 

3 ICREA-Instituci6 Catalana de Recerca i Estudis Avangats, 08010 Barcelona, Spain 



Abstract. 

We employ the exact diagonalization method to analyze the possibility of generating 
strongly correlated states in two-dimensional clouds of ultracold bosonic atoms which 
are subjected to a geometric gauge field created by coupling two internal atomic states 
to a laser beam. Tuning the gauge field strength, the system undergoes stepwise 
transitions between different ground states, which we describe by analytical trial wave 
functions, amongst them the Pfaffian, the Laughlin, and a Laughlin quasiparticlc 
many-body state. The adiabatic following of the center of mass movement by the 
lowest energy dressed internal state, is lost by the mixing of the second internal 
state. This mixture can be controlled by the intensity of the laser field. The non- 
adiabaticity is inherent to the considered setup, and is shown to play the role of 
circular asymmetry. We study its influence on the properties of the ground state 
of the system. Its main effect is to reduce the overlap of the numerical solutions with 
the analytical trial expressions by occupying states with higher angular momentum. 
Thus, we propose generalized wave functions arising from the Laughlin and Pfaffian 
wave function by including components, where extra Jastrow factors appear, while 
preserving important features of these states. We analyze quasihole excitations over 
the Laughlin and generalized Laughlin states, and show that they possess effective 
fractional charge and obey anyonic statistics. Finally, we study the energy gap over 
the Laughlin state as the number of particles is increased keeping the chemical potential 
fixed. The gap is found to decrease as the number of particles is increased, indicating 
that the observability of the Laughlin state is restricted to a small number of particles. 
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1. Introduction 

A spectacular progress on the manipulation and control of cold atomic clouds has been 
achieved since the first experimental realization of a Bose-Einstein condensed system 
by Anderson et al. jl] and Davis et al. [2] nearly simultaneously in 1995. On the one 
hand, these systems provide us with a toolbox for studying the principles of quantum 
mechanics, as has for instance been done in experiments involving the interference of 
two condensates leading to the manifestation of coherence [3]. On the other hand, 
cold atoms have been used in order to simulate interesting phenomena appearing in 
condensed- matter physics jl], like the Mott-insulator to superfluid transition [5]. 

A severe restriction, however, stems from the atom's electro-neutrality, which 
hinders a direct implementation of phenomena involving electromagnetic forces. An 
important example of the latter is the physics of the integer and fractional quantum 
Hall effect, occurring in two-dimensional interacting systems under the presence of a 
strong perpendicular magnetic field. The ground states of fractional quantum Hall 
(FQHE) systems are highly correlated states, like the Pfaffian state proposed by Moore 
and Read [6], or the celebrated Laughlin state [7], whose bosonic analogs are found to be 
the exact eigenstates of a Hamiltonian with a three-body (3B), or two-body (2B) [El Ej 
contact interaction, respectively. In addition, the quasihole excitations over the Laughlin 
state have fractional effective charge and fractional statistics [in]. Excitations of the 
Pfaffian state may even obey non-Abelian braiding statistics [HI [12] , which makes them 
interesting from both fundamental and technological points of view [13j . 

Several proposals of experimental routes to obtain these types of states have 
appeared in the literature, ranging from the use of rotating traps to simulate magnetic 
fields acting on charges to the use of laser-beam configurations acting on atoms with 
several internal states HH HH1 UHl HZl HH HH1 [201 [21]. In this article we analyze 
the appropriate conditions to realize some strongly correlated states within a simple 
configuration of a single laser beam shining on a cloud of cold atoms with two internal 
states. If the internal dynamics of the atoms, governed by the Rabi frequency of 
the atom-laser coupling, is fast enough with respect to the slow variation of center 
of mass position, then the atoms evolve adiabatically. They remain in one definite 
space-dependent superposition of the internal bare states, and the accumulation of Berry 
phase [22] during its movement mimics an effective magnetic field [231 121] • An important 
goal is to go further and analyze the effect of a slight amount of non-adiabaticity, which 
we treat in a perturbative way. 

Through the controlled variation of external parameters, different strongly 
correlated states appear in the spectrum, i.e. a Laughlin-like state, a Pfaffian-like state 
and the quasiparticle-like state obtained from the Laughlin. Our main aim is to map the 
regions in parameter space where the exact, numerically obtained, ground state (GS) has 
a large overlap with explicit analytical expressions provided for these relevant strongly 
correlated states. To remain close to possible experimental implementations, we study 
the effect of small perturbations that create non-adiabaticities which are unavoidable 
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for finite values of the laser intensity. The small perturbation produces a deformation 
of the atomic cloud preserving most of the notable properties of the original states, e.g. 
entropy, internal energy, or anyonic character of excitations. These slightly asymmetric 
ground states are well represented by generalized analytic wave functions. 

Focusing additional laser beams on the atomic cloud, it is possible to pierce holes 
into the system. As long as the asymmetric perturbation is small, the resulting states 
can be well described by an analytical quasihole wave function, which can be obtained 
from the Laughlin or the generalized Laughlin wave function. In both cases, the effective 
charge and statistical phase of the quasihole excitations are found to attain fractional 
values, demonstrating the possibility of observing anyons within the proposed setup. 

In this work we concentrate on the physics of few-body systems, independently 
of their attainability in the thermodynamic limit, which is beyond the scope of the 
present article. This approach is meaningful as there are nowadays a number of 
groups able of dealing with small bosonic clouds using several techniques [251 126] . 
In particular, Ref. [26] has presented experimental evidence for the production of 
quantum states of fractional quantum Hall type for small atom systems (N<10). These 
experimental developments have triggered a number of theoretical proposals focusing on 
the production of strongly correlated quantum states in small atomic clouds [161 EE ISHj • 
We analyze the dependence on the number of particles of some of our main results. First, 
we show that observing the anyonic character of quasihole excitations becomes possible 
when N > 5. Second, we show that at fixed chemical potential the bulk gap of the 
Laughlin state is decreased as N is increased, indicating that few-body systems provide 
the best scenario for producing the bosonic Laughlin state. 

The paper is organized as follows. In Section [2] we present the system and 
derive an effective Hamiltonian to describe it. In Section [3] we provide the analytical 
expressions of the relevant fractional quantum Hall states and their generalizations. Our 
results are presented in Section H] and Section [5] in the adiabatic/symmetric and non- 
adiabatic/asymmetric cases, respectively. In Section [6] we study the fractional charge 
and anyonic statistics of quasiholes, produced by means of additional lasers. In Section [7] 
we analyze the behavior of the energy gap above the Laughlin state as N increases, and 
its evolution as a function of the system's parameters. Finally, in Section [HI we present 
our conclusions. 

2. Description of the system 

We consider a setup to produce artificial gauge fields in a small cloud of ultracold atoms 
closely following the configuration described in Ref. [21]. The system is confined by 
harmonic traps. The confinement in the z direction is assumed to be strong enough to 
achieve effectively a two dimensional system. The cloud is illuminated by a single laser 
beam with wave number k and frequency cul, which propagates in the y-direction and 
is close to the resonance with a transition between two internal atomic states, ul —^a- 
The interaction between the electric field of the laser and the induced electric dipole is 
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Hs » = 2M + V & + ~2 { e^sinfl -cos J ' (3) 



modeled by the atom-laser Hamiltonian, which in the rotating-wave approximation and 
in the rotating frame is given by [27J 12H] , 

H AL = E s |g> (g| + (E e - hu L ) |e) (e| + |e> (g| + |g) (e| (1) 

where E s and E c are the energies of the bare atomic ground and excited state and 
Qq is the Rabi frequency, which is proportional to the laser intensity No spontaneous 
emission of photons from the excited state is considered. This assumption is justified 
as long as the lifetime of the excited state is longer than the typical duration of an 
experiment. Lifetimes of several seconds, as found for Ytterbium or some Alkaline earth 
metals, should be sufficient. 

In order to obtain a non-trivial gauge potential from the coupling in Eq. ([1]) , we still 
have to make it dependent on x. This can be achieved via a real magnetic field, which 
by Zeeman effect makes the internal energy levels vary linearly with x. Introducing a 
parameter w, setting the length scale of this shift, we have, 

E, = -^*, E, = Hv A+ ^*. (2) 
2 w 2 w 

Then, the single particle Hamiltonian is given by 

p 2 Ml f cos 6 e^sinO 

2A 

where the third term is the atom-laser Hamiltonian represented in the {|e) , |g)} 
basis. Here, M is the atomic mass, f2 = Sl \/ 1 + x 2 /w 2 , sin 6 = w/yfw 2 + x 2 , 
cos# = x/Vw 2 + x 2 , and = ky. V(r) is the trap potential that will be fixed below. 

Up to this point, the Hamiltonian is given by physically measurable parameters. 
In the next step, we choose particular expressions for the single particle states that 
diagonalize the atom-laser Hamiltonian. This corresponds in fact to the selection of 
a specific gauge for the vector and scalar potentials that will drive the center of mass 
dynamics, as will be shown below. We consider the eigenf unctions of Hal given by, 

/ n p i<t>/2 \ / _ o i0/2 \ 

in the {|e) , |g)} basis, where C = cos 0/2, S = sin#/2, and G = The atomic state 
can be then expressed as, 

x (f) = ai(f) ® + 02(f) (8) ^2) (5) 

where Oj captures the external dynamics and provides internal degree of freedom. 
The single-particle Hamiltonian is then expressed in the {|^i) , ^2)} basis as 

\n 2 i tiii j 

acting on the spinor [ai(r), a-i(ff\- Defining a vector potential A, 

y x x 



A(f) = hk 



Aw' Aw 2\f. 



x L + w 



(7) 
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and a scalar potential U, 



h 2 w 2 



U{r) 



8M (x 2 + w 2 ) 



(8) 



we obtain, 



p — A(r) 



and 



H. 



2M 



p + A(r) 



22 



+ U(r) + V(r) + 



+ U(f) + V(f) 



(9) 



(10) 



2M v ' v ' 2 

which are the Hamiltonians driving the external dynamics of atoms being in the internal 
state and ^2), respectively. By expanding the terms up to second order in x and 
y, which is justified by choosing w to be larger than the extension of the cloud, we find 
that the energy distance between these two manifolds is given by HQq. For convenience, 
we make the Hamiltonian element for the low energy manifold, if 22, independent of 
VIq by adding the constant term to the diagonal of H sp . Further we note that 
with the explicit selection of Eq. (j4]) and for x, y -C w, A is in the symmetric gauge: 
^ ~ T~{Vi —x)- This allows for making i?22 fully symmetric by a proper choice of the 
trapping frequency Eq. (flO]) then takes the form 

.2 *-£m 



H- 



V 



22 



and we can take the advantage from the knowledge of its single particle eigenfunctions, 
namely the Fock-Darwin states [32]. This final expression is formally equal to the 
Hamiltonian driving a system of charges trapped by a harmonic potential of frequency 
u±, under a constant magnetic field along the ^-direction, or equivalently, a system of 
neutral atoms trapped by a rotating harmonic potential of frequency u±, expressed in 
the rotating frame of reference [HU [34] . 

While the equivalence to the rotating case holds for H22, we stress that it does not 
for the system's behavior in the upper manifold described by Hn. However, due to the 
off-diagonal terms in H sp , both manifolds are coupled. Typical expected values of if 12 
and H 2 i are of the order of the recoil energy Er = which gives the scale for the 
kinetic energy of the atomic center-of-mass motion when it absorbs or emits a single 
photon. If we consider MIq ^> Er, this coupling is small, and we can restrict ourselves 
to the low energy manifold. Namely, we are in the situation where the internal dynamics 
is much faster than the center of mass motion and can follow the external variations in 
a quasiadiabatic way [29] . 

To go beyond the adiabatic approximation, we consider the influence of the high- 
energy manifold as a small perturbation. Using the procedure appropriate to systems 
that show two significantly different energy scales as explained in Ref. [27], and detailed 
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Appendix A, we calculate an effective Hamiltonian up to second order in the 



perturbation, which reads 



H ™ ~ H22 ~ ~W (12) 

where the explicit expression for the perturbative term ^21-^12/^0 is given in 



Appendix A Though mathematically more involved and physically richer, this term 



is reminiscent of the anisotropic potential that is applied to set an atomic cloud 
in rotation. Usually, the expression used to model the stirring laser is given by 
aMu±(x 2 —y 2 ) [301 El] , where a measures the strength of the deformation. Similarly the 
term f/21-^12/^0 can only produce changes of L of AL = 0, ±2, ±4 (see Appendix A ). 
In what follows we will identify 'deformation' with 'coupling' indicating that a large 
connection with Hn implies deformation since Hu is not cylindrically symmetric. 

The many-body Hamiltonian which we finally will deal with is given by adding a 
contact interaction term to the effective Hamiltonian from Eq. ( [121) : 

N h 2 

H = J2 H%{i) + ^ E ^ - ? i) ( 13 ) 

i=l i<j 

where g is a dimensionless parameter fixing the contact interaction strength. From now 
on we will consider gN = 6 in the numerical calculations. The Hamiltonian H acts only 
on the low energy manifold, which effetively is perturbed by the second manifold. The 
important two parameters in H are the dimensionless ratio 77 = 4A /^ and the degree 
of the perturbation given by Qq. The expression rjui± plays the role of the rotating 
frequency, and also the effective magnetic field strength B is proportional to 77: 
hk 

Bo = ^ = 2Mui»). (14) 

The largest possible value of 77 is given by 1 in order to keep the system confined. 
We consider that the artificial magnetic field is strong enough to work exclusively in 
the lowest Landau level regime. To this end, our parameters must fulfill the following 
conditions: the energy difference between Landau levels is larger than both, the kinetic 
energy of a single particle within a Landau level and the interaction energy per particle. 
Note that the strength of the atom-laser coupling, characterized by Qq, is different from 
the strength of the magnetic field, characterized by rj. Because the magnetic field has a 
geometric origin, 77 is independent of the atom-laser coupling, as long as the adiabatic 
approximation is justified. 



3. Analytical many-body states 



In this section we give an overview of the analytical wave functions discussed in the 
literature which will turn out to be relevant for describing our system. 
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3.1. Laughlin state 

The well-known Laughlin state has the analytical form [7j (35J [36] , 

9c( Zl , . . . , z N ) = Af c H( Zi - % .)i/V^/»i , (15) 



i<j 



where Mc is a normalization constant, z = x + iy and Aj_ = ^/h/Muj±. The inverse 
of the exponent, v, fixes both the density of the system and the symmetry of the wave 
function. For bosons, \ jv must be even. The Laughlin regime discussed here will be at 
half filling, so for the rest of this paper we set v — -x. 

The analysis of the squared overlap | (\I/£|GS) | 2 of the Laughlin state with the exact 
ground state as a function of the artificial magnetic field strength rj and the atom- 
laser coupling Q shows that the overlap gets reduced as hfl /E R decreases, even for 
large values of rj. Larger overlap with the ground state can be obtained by adding 
an admixture of the Laughlin state with additional Jastrow factors that allow for an 
increase of total angular momentum, which in the Laughlin state is given by the integer 
L = N ^~^ ■ Based on these observations, in Ref. [2T] an analytical ansatz for the 
ground state in the Laughlin-like region [j] was proposed, 

^gc = a^c + PVn + 7^£2 , (16) 
from here on referred to as the generalized Laughlin (GL) state, with \&cx = A/1 ^ c ■ 
Y^=izh ^C2 = M 2 {^C2 - and § £2 = Af 2 ^c ■ Y^<jZi z V such that 

we ensure (^c\^a) — and (^al^Cj) — This ansatz involves components of 
angular momentum L = N(N — 1) and L = N(N — 1) + 2, and yields zero contact 
interaction energy. The values of a, (3 and 7 are computed as a = (\& £|GS) / VXf \ 
= (* £1 \GS)/Vtf, and 7 = (Vc2\GS)/VJ7, with J\f = |(^|GS)| 2 + |(^ £ i|GS)| 2 + 
|(^£ 2 |GS)| 2 . 

3.2. Pfaffian (Moore-Read) state 

While the Laughlin and the generalized Laughlin state turn out to be good trial states for 
strong magnetic fields rj < 1, for smaller field strengths the Laughlin quasiparticle (LQP) 
state and the Pfaffian (Pf ) state become relevant. The Pfaffian state has L = N(N—2)/2 
for even N, and L = (N — l) 2 /2 for odd N and its analytical expression reads, 

9 v =tf vS Pf([z])J[(z i -z j ) (17) 

i<j 

with jVpf a normalization constant, and 



\Z\ — Z 2 ) (z 3 — Z4) (zn-i — Zn)_ 

where A is an antisymmetrizer of the product. As explained in Ref. [HI E], the Pfaffian 
state can also be computed as, 

V P = S Y[ {Zi - z 3 f J] (z k - Zl ) 2 (19) 



i<j£ai k<l£ci2 

I Defined as the parameter domain where (L) > N(N — 1). 
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where <j\ and a 2 are two subsets containing N/2 particles each if N is even, and (A r +l)/2 
and (N — l)/2 if N is odd. S symmetrizes the expression. The Pfaffian state has been 
shown to be the lowest energy eigenstate of a Hamiltonian which contains only three- 
body contact interactions [6| [20]. 

H int = S(zj - Zj)S(zj - Zk) , (20) 

i<j<k 

and, similarly to the Laughlin with the two-body interaction, it has zero three-body 
interaction energy. Remarkably, however, as shown in Refs. [El EE] for a two-body 
interaction Hamiltonian and for some values of rj there is a sizeable overlap between the 
ground state of the system and this state. 

As we did before for the Laughlin state in Eq. f lT6|) . we can define a generalized 
Pfaffian (GPf) state as 

tyg V = a^-p + P^ V1 + 7^p 2 , (21) 

with = J\f P1 ^ P ■ Yn=i z h ^V2 = N V2 (^>v2 ~ (*pi|*7>2)*7>i), and ^ V2 = 

Mv2^v ■ J2i<j z i z j- Again, the parameters a, (3, and 7 are fixed to maximize the 
overlap of the numerical ground state with ^g-p. 

3.3. Laughlin- quasiparticle state 

The Laughlin-quasiparticle state arises from the Laughlin state by increasing the density 
at the origin, decreasing its angular momentum, L qp = N(N—1)—N. The latter formula 
holds if the quasiparticle is at the origin. Otherwise, it also carries angular momentum 
and the total expected value of the angular momentum of the system is no longer an 
integer. The wave function is written as, 

^ qp = AU£, O (d Zl - • • ■ (d ZN - 0^£ , (22) 

with A/* q p(£, £*) a normalization constant that depends on the position £ and £* of the 
excitation. Also for the Laughlin-quasiparticle state we define a generalized version 
(GLQP), built up from the same Jastrow factors used in Eq. fll6j) . i.e. 

VgCw = a*£qp + /3*£qpl + 7*£q P 2 , (23) 

With Vcqpl = A/" £qpl ^ qp • ^=l Z h #£qp 2 = A^ qp2 (^ £qp2 - (^£qpl | ^£qp 2 > ^£qpl) , and 
*qp2 = Afc qP 2 ^£qp " Z t Zj . 

3.4- Laughlin- quasihole state 

Alternatively to increasing the homogeneous density of the Laughlin state locally, one 
might also decrease it by piercing a hole in the atomic cloud. Formally this is achieved 
by introducing an additional zero into the wave function, multiplying it with n?(£ — z i)- 
The resulting quasihole state pierced at £ reads: 

N 

*£qh = A^ qh (£, C) " > ( 24 ) 

1=1 
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where A/£ q h (£,£*) is a normalization constant, which explicitly depends on the position 
of the quasihole. To test the anyonic nature of the quasiholes, we need to do the same 
operation twice, such that the presence of two quasiholes is described by, 

N 

^£ 2qh = A^ 2qh (6, Ci, 6, £) II^ 1 - z ^ - • ( 25 ) 

i=l 

The state with one quasihole has a fixed total angular momentum, which is N quanta 
above the ground state, if the quasihole is at the origin, £ = 0. For off-centered quasihole 
positions, the average angular momentum is slightly reduced and non-integer. The state 
with two quasiholes has an angular momentum close to 2N quanta above the ground 
state. 

We may also apply the same operation to the generalized Laughlin state and define, 

N 

*g£ qh = A/" g£qh (£, t) - *) + + ( 26 ) 

i=l 

for the state with one quasihole, and 

^gC2qh = A/g£2qh(£l> £l, £2, £2) 

TV 

x - - **) + + ( 2? ) 

i=l 

for the state with two quasiholes. As explained in Sec. El we always find a, (3 3> 7, thus 
in practice we will consider always 7 = 0. 

4. Results for the adiabatic/symmetric case, ifff = H22 

In the symmetric case, in which the perturbation if 21-^12/ (^0) is not included, and for 
N = 4, four distinct regions are detected depending on the value of 77 as previously 
obtained in Ref. [H [16]: Condensed, Pfaffian, Laughlin-quasiparticle and Laughlin 
regime. We analyze them in the following Figs. [UandEJ 

The first region corresponds to a fully condensed system with zero angular 
momentum and vanishing one-body entanglement entropy |§|, see Fig. [TJ The ground 
state can be well described by a wave function given by — •N"c er ^ Z V 2X ±- being Ac 
a normalization constant. 

Figs. Q] (middle panel) and |2] show that the first excitation is of quadrupolar 
character up to rj = 0.7. Lower values of 77 lie beyond our lowest Landau level 
approximation, where a larger Hilbert space including more Landau levels has to be 
considered. 

At the critical value 771 = 1 — gN/(8ir) ~ 0.76, (gA^ = 6), a degeneracy between 
states with L — 0, 2, 3, and 4 occurs, see Fig. [2j At this r\\ a state with broken 

§ The entropy is defined here from the one-body density matrix, as S = — ^rii ln(rii), where nj are 
the eigenvalues of the one-body density matrix. For a more detailed discussion of the entropy we refer 
for instance to pH] , 
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Figure 1. (Color online) (upper panel) Interaction energy in units of hui± (black 
circles) and one-body entanglement entropy (red squares) of the ground state as a 
function of r/. (middle panel) Angular momentum in units of h of the ground state and 
of the first excited state as a function of 77. (lower panel) Squared overlap between the 
ground state of the system and the exact Laughlin, Pfaffian and Laughlin-quasiparticlc 
states. The plots corresponds to the case £Tff = i?22- 



symmetry, combination of several L-eigenstates, is precursor of the nucleation of the first 
vortex state. For increasing 77 the ground state recovers the cylindrical symmetry and 
the angular momentum is L = 4. All this phenomenology can be inferred from the Yrast 
line displayed in Fig. [3j The Yrast line is constructed by plotting the interaction energy 
contribution of the lowest energy state for each L. From this line, the addition of the 
kinetic energy, which reads (up to a term independent of L and rj) E kin = (1 — rj)Lfkj±, 
produces the total energy with its minimum at the angular momentum of the ground 
state, Lqs, as exemplified for rj = 0.85 and 0.94 in Fig. [3j This is a general behavior for 
any N. 

As rj is increased above rj\, the ground state with L = 4 becomes correlated as 
marked by the sizeable entanglement entropy (see Fig. [1]), meaning that it is not fully 
condensed. Remarkably, its squared overlap with the Pfaffian state is much larger 
(about 0.9, see Fig. H]) as already noticed in Ref. [H [16], than its squared overlap with 
the one- vortex-state (about 0.47), given by = A/ivx YliLi Zie~^ iZ ^ 2X± (being A/i vx a 
normalization constant). If the ground state in this region is modeled by the eigenstate 
of the one-body density matrix with the highest occupation which plays the role of the 
order parameter in a mean field approach, it shows a well defined vortex at the center 
of the density distribution. The phase of the order parameter changes by 2tt, indicating 
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Figure 2. (Color online) Energy difference, in units of Hui±, between the ground state 
and its first excitation as a function of rj. The large blue numbers correspond to the 
value of L for the ground state. The small numbers quote the value of L of the first 
excited state. TJff = H^i- 



the existence of a vortex. However, this function is a poor representation of the ground 
state since the non-condensed fraction is significant. This region has three different kind 
of excitations, L = 3, 6 and 8 as can be seen in Fig. |2j The latter has a large overlap 
with the Laughlin-quasiparticle state as can be seen in Fig. HI 

For 0.92 < T] < 0.96, the ground state has L = 8, a higher entanglement entropy, 
and a smaller interaction energy. The ground state has a large overlap with the Laughlin- 
quasiparticle state, as shown in Figs. [Hand HI This quasiparticle region has four different 
excitations, L = 4, 6, 9 and 12. The L = 9 corresponds to a center of mass excitation of 
the system as dictated by Kohn's theorem [321 [37] . 

Finally, for 77 > 772 — 0.96, the ground state wave function is the Laughlin wave 
function, with zero interaction energy. Its excitations are the Laughlin-quasiparticle 
L = 8 and an a center of mass excitation, Kohn mode, with L = 13, whose analytical 
form is, ^ = N(z 1 + z 2 + z 3 + z 4 )^ c ■ 

5. Effects of the non-adiabaticity /asymmetry, iJ^f = -^22 — -£^21-^12/(^0)- 



As discussed in Section [2] and Appendix A, the considered setup can be mapped onto 



a symmetric Hamiltonian, H 2 2, equivalent to the one of rotating atomic clouds in 
symmetric traps, plus a term, ^21-^12/(^0) whose importance can be controlled by 
tuning the laser coupling, Q . As discussed in Ref. [21], the first effect of the perturbation 
in the Laughlin-like region is to increase the angular momentum of the ground state by 
populating the states ^ci, defined in Eq. ( TIB"]) . One can consider fairly small coupling 
MIq/Er ~ 40 and still get Laughlin-like ground states of the form of Eq. (fTBT) . which 
retain most of the known properties of Laughlin states, namely, a large entanglement 
entropy and vanishing interaction energy [21]. Now we extend the previous study to the 
effect of the perturbation on the Pfaffian and Laughlin-quasiparticle regions. 

In Fig. [5] we show the squared overlap between the ground state and the three 
original correlated states (left panel) and their generalized versions (right panel) 
identified as the generalized Pfaffian, generalized Laughlin-quasiparticle (GLQP) and 
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Figure 3. (Color online) Yrast line for N — 4, solid red circles, which corresponds 
to the interaction energy contribution of the lowest eigenstates for each value of L. 
The triangles and diamonds depict the sum of the interaction energy and the kinetic 
contribution for i] = 0.85 and i] = 0.94, respectively. The arrows mark the value of L 
which corresponds to the GS in each case. The energies are given in units of thvj_. 
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Figure 4. (Color online) Squared overlap between the Pfaffian, the Laughlin- 
quasiparticle and the Laughlin states and the ground state (upper panel) and the 
first excited state (lower panel). The condensed region has been omitted. = i?22- 



generalized Laughlin (GL), see Eqs. ffl~6]) . (|2ip . and (|23|) . It turns out that in all 
three cases the state which is proportional to 7 is much less populated than the 
states proportional to a and (3. We thus neglect the contribution of this term, for 
simplicity. As shown in Fig. |5l overviewing all three regimes, the largest improvement 
by using the generalized versions occurs in the Laughlin region. Here, the total angular 
momentum increases continously with 7] [21], leading to substantial occupation of the 
state proportional to j3. 
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Figure 5. (Color online) (left panel) Squared overlap between the ground state and 
the original strongly correlated states considered, namely, the Pfaffian, Laughlin and 
Laughlin-quasiparticle states as a function of r] for hflo/E^ = 40 and 100. (right panel) 
Squared overlap between the GS and the generalized correlated states considered, GPf, 
GL and GLQP as a function of 77 for MIq/E^ = 40 and 100. 
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Figure 6. (Color online) Squared overlap between the ground state and the Pfaffian 
and generalized Pfaffian states defined in the text, upper and lower panels, respectively, 
as a function of r\. The different lines correspond to different values of HUo/Eb,. The 
solid line is obtained with iJ^f = #22 • 



Detailed information about the effect of the perturbation is shown in Figs. \5f7\ and 
[8] for each of the three regions seperately. First, in Fig. [6] we consider the overlap with 
the Pfaffian and GPf, exploring fairly low values of KIq/Er. Lower values of MIq/Er 
require the consideration of higher order terms in the expansion oftf 2 f not included 111 
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Figure 7. (Color online) Squared overlap between the ground state and the Laughlin- 
quasiparticle (LQP) and generalized Laughlin-quasiparticle (GLQP) states denned in 
the text, upper and lower panels, respectively, as a function of r\. The different lines 
correspond to different values of KIq/Er. The solid line is obtained with ifff = -^22- 



our calculations. There one can see how by decreasing the value of WIq, the rj which 
provides maximum overlap becomes smaller. Thus, while in the symmetric case, the 
only region with non-negligible squared overlap with the Pfaffian was 0.75 < r\ < 0.92, 
with, e.g hflo/En = 20, the region is roughly displaced to 0.73 < rj < 0.89 

Also, the squared overlap with the Pfaffian gets reduced, going from around 0.9 in 
the symmetric case, to 0.7 for KIq/Er = 10. As occurred with the Laughlin, the main 
effect of the perturbation is to populate states which are of the GPf type, that is, a 
Pfaffian core with appropriate Jastrow factors. The GPf state has a large overlap with 
the ground state, of the same order as the Pfaffian itself with the symmetric ground 
state. Interestingly, large values, > 0.8, of the squared overlap with the GPf state can 
be found already for HVLq/Er > 20, which is relevant from the experimental point of 
view as it increases the window of observability. 

A similar behavior is found when studying the squared overlap of the Laughlin- 
quasiparticle state with the exact ground state of the system. As shown in Fig. [7] the 
region with sizeable overlap with the Laughlin-quasiparticle state gets shifted towards 
lower values of 77 as we decrease f2 , peaking at 77 = 0.85 for T&Iq/Er = 10. Also, a 
sizeable overlap is found with the GLQP state. It is however clear, that large values of 
the squared overlap, >0.8, can only be found for values of KVLq/Er > 30. 

In Fig. |8] we present the corresponding figure for the case of the Laughlin and 
generalized Laughlin states. First, we note that again the region where the L and 
GL are most populated gets shifted towards lower values of rj as we decrease the 
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value of KVLq/Er. For instance, the maximum value is obtained around rj = 0.91 for 
MIq/Er = 15, and this maximum is of 0.4 in the case of GL. Squared overlaps larger 
than 0.8 are only obtained for MIq/Er > 40. Squared overlaps larger than 0.5 can 
however be obtained with HVLq/Er as low as 20 for rj ~ 0.92. 
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Figure 8. (Color online) Squared overlap between the ground state and the Laughlin 
(L) and generalized Laughlin (GL) states defined in the text, upper and lower panels, 
respectively, as a function of r\. The different lines correspond to different values of 
hQ /E R . The solid line is obtained with iJ|f = H 2 2- 
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6. Fractional charge and anyonic statistics of quasihole excitations in the 
Laughlin regime 

An important property of the strongly correlated iV-body states discussed in the 
previous sections are their excitations which might behave as particles with fractional 
charge and obey anyonic statistics, as is the case for the quasihole excitations over the 
Laughlin ground state [lOj 115] . 

An experimentally feasible way of creating quasihole excitations in our system is 
by focusing a laser beam onto the atomic cloud. This can be described by adding the 
following potential to the single-particle Hamiltonian of Eq. ( 1T3|) |15j : 

i 

where / is the laser intensity, and £ is the position onto which the beam is focused. With 
two such potentials we should be able to create states with two quasiholes, according 
to Eq. (|25|) . In the following, we will first study quasiholes in the Laughlin state, which 
can be created when the system is in the adiabatic regime. Then we will also consider 
a slightly non-adiabatic situation, where we find quasiholes in the GL state, as defined 
by Eq. fl26|) . We analyze the quasiholes in both the Laughlin and the GL state with 
respect to their fractional character. 



6.1. Quasiholes 1 properties in the adiabatic case 

As already discussed in Section HJ for h~Q 3> E R , the system's ground state squared 
overlap with the Laughlin state is effectively one, above the critical field strength 772. 
Now we consider the system with the additional term fl28|) and find that there is also 
a region of 77 where the overlap of the ground state of the system and the analytical 
quasihole state is effectively 1, see Fig. [9j This shows that the potential of Eq. (|28|) 




0.99 



Figure 9. (Color online) Squared overlap of exact ground states (N = 4) with 
Laughlin state (solid) , Laughlin state with one quasihole (red) , and Laughlin state with 
two quasiholes (blue). The quasiholes are created by a laser with intensity / = 

(dotted lines) and / = 30^- (dashed lines). 
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Figure 10. (Color online) Moving one or two quasiholes at fixed radial positions on 
circles around the origin allows for determining the effective charge and the statistical 
phase angle. The quasiholes are always fixed at a radial position R = 1 (in units of 
Ax). 



is able to produce quasiholes described by ( )24|) . Similarly, adding two such lasers we 
also find a region of r] where the overlap between the ground state and the analytical 
state with two holes, Eq. ( )25|) . is very close to 1. However, we notice that the values 
of r) at which the overlap for one or two quasiholes reaches 1 differ from each other, 
and both are found for values larger than r] 2 , see Fig. [91 These features are essentially 
independent of the laser strength /, for sufficiently large /. 

The most interesting property of these excitations is their fractionality, i.e. 
fractional charge and statistics. To study the fractional charge of the quasiholes, we 
first note that in our electro-neutral system subjected to an artificial magnetic field, 
there exists the analogue of an electric charge which can be defined via the behavior 
of a particle or quasiparticle evolving within the gauge field. To this end, we consider 
the phase a quasihole picks up while being adiabatically displaced following a closed 
trajectory. The general expression for the Berry phase on a closed loop C is given 
by [3H] 



7c 



(29) 



with \^(R)) the state of the system, characterized by a parameter R, which in our case 
is the position of the quasihole. For simplicity, we now assume that the quasihole is fixed 
at a radial position |£| = R\±, but is moved along a circle centered at the origin, see 
Fig. HO, parameterized by the angle <fi. This is sufficient to test the fractional behavior. 

For general contours, one can extract the acquired phases from the normalization 
factor of the quasihole state, as described in [39]. For our circularly symmetric contour, 
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however, the situation is simpler, as we can re-write Eq. ( 12 9 p as 




(30) 




Here we note that, due to the circular symmetry of the Laughlin state, the integrand 
does not depend on the angular position of the quasihole. The function f(R) can be 
calculated by decomposing the Laughlin quasihole state into Fock-Darwin basis [40J, 
which we have done analytically for particle numbers up to 6. For compactness, we 
explicitly give here only the result for N = 4: 



If we assume that the quasihole is moved sufficiently close to the center, i.e. R < 1, we 
can expand this expression in R and find, 7 C = 2tt(0.473426 j R 2 +0.0242202 j R 4 +C( j R 6 )) w 
ttR 2 . Thus, the acquired phase is approximately given by the enclosed area in units of 
Ai. ^ 

To obtain the effective charge of the quasihole, we must compare this result with 
the geometric phase acquired by a particle moved along the same closed contour due 
to the gauge field. In the Laughlin regime, where r\ ~ 1, we find with Eq. (1141) 
Bq pa 2^-, thus the acquired phase (p is two times the enclosed area in units of Aj_, 
i.e. tp = j i Bq(RX±) 2 tt pa 2R 2 n. From this follows the effective charge of the quasihole 
to be g e ff = — ~ 0.47, close to the expected value for the Laughlin state at half filling 
in the thermodynamic limit, 1/2 [10]. We have performed a similar study for A^ = 5 
and N = 6, finding that for N = 5, the effective charge is about 0.48, and for N = 6 it 
is found to be 0.49, i.e. by increasing the particle number the value 1/2 is approached. 

6.1.1. Fractional statistics 

To prove the fractional statistics of the quasihole excitations we now consider the 
system with two quasiholes at £i = l^ile** 1 and £2 = l&le , which we assume to sit on 
opposite radial positions, i.e. |£i| = |^| = -RAj_ and <p2 — 4>i = 71 ■ We now consider the 
simultaneous adiabatic movement of the two quasiholes on two half circles, in such a 
way that, at the end, the quasiholes interchange position (see Fig. [T0|) . This differs from 
a more common setup to test the statistical angle, where one quasihole is fixed in the 
center, while the other is encircling it, but it has the advantage that it maximizes the 
distance between the two quasiholes. Note that in Fig. [TUJ the radial position is chosen 
at R — 1, i.e. the distance between the center of the quasiholes is 2Aj_, which seems to 
be the minimum distance needed for not having a significant overlap (< 10%) between 
the two quasiholes, see Fig. HIT a). On the other hand, in this small system of just four 
particles, larger radial positions lead to quasiholes overlapping with the systems edge. 

The total phase picked up during the described movement should be the sum of 
the phase picked up by one quasihole moved along a circle plus a phase factor due 
to the interchange of the two quasiholes. Again the phase gradient turns out to be 



7c = 2nf(R) = 2tt 



4 (10128i? 2 + 5313# 4 + 1659i? 6 + 553i? 8 ) 



(31) 



85572 + 40512/2 2 + 10626i? 4 + 2212/? 6 + 553i? 8 
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Figure 11. (Color online) (a) Squared overlap between the quasihole wave functions, 
|(^ , £qh(Ci)| v I'£qh(C2))| 2 at opposite angular positions, <\>\ — <f>2 — 7T, as a function their 
distance to the center, |£i| = |^| = i?Aj_. (b) Statistical angle of two quasiholes at 
opposite angular positions and radial position RX± . In both panels we present results 
for N = 4 (black solid), N = 5 (red dotted), and N = 6 (green dashed). 



independent from the angular position, but is described by a different function / of the 

radial coordinate: 

ft p>\ - 8 (286812CLR 4 + 461616i? 8 + 25242i? 12 + 553i? 16 ) 

^ ' ~ 41660640 + 11472480/? 4 + 923232i? 8 + 33656.R 12 + 553i? 16 ' ^ ' 

The statistical phase angle is thus, 

¥>(i*W = / /(^)# - / f(R)d<f> = 2nf(R) - nf(R) . (33) 
Jo Jo 

It is shown, as a function of R, in Fig. [11] (b). First, as expected the statistical phase is 
zero if both quasiholes are in the same position, and it increases linearly as the distance 
between the quasiholes is increased. This linear behavior then saturates once the overlap 
between the two quasiholes, |(^£ q h(^i)|^'£qh(^2))| 2 , drops below 0.1, and remains mostly 
constant around ir/2. By increasing the number of particles N, the phase angle become 
less dependent on R once the two quasiholes do not overlap, meaning that it becomes 
a robust property of the quasiholes. It stabilizes around the expected value of tt/2. At 
larger distances R, the system's edge starts to play a role. 



6.2. Non-adiabatic effects on the properties of quasiholes 

To study the fractionality of quasihole excitations in the non-adiabatic case we will again 
profit from the generalized analytical representations used to describe the ground state 
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Figure 12. (Color online) Left: Squared overlap of exact ground states (N = 4) with 
generalized Laughlin state (solid black) , generalized Laughlin state with one quasiholc 
(red), and generalized Laughlin state with two quasiholes (blue). The quasiholes are 
created by a laser with an intensity / = 30-^^ (dotted lines) and / = 50-^- (dashed 
lines) at £i = X± and £2 = — Aj_. The Rabi frequency is M^o = 100-Er. Right: The 
values of the variational parameter (3 for a given r\ are similar for all three states. 



of the system. Following the discussion in the previous section we compute the squared 
overlap of the ground state obtained with no-, one-, and two-extra lasers piercing holes 
into the system. First, we find a significant squared overlap for the slightly perturbed 
case at HQq = IOOEr, see Fig. [12j As occurred in the adiabatic case, a large overlap with 
the analytical one- and two-quasihole states appears only at higher field strengths than 
the one at which the generalized Laughlin state is reached. Our study of the properties 
of quasiholes in the non-adiabatic case will be restricted to the parameter domain where 
a fair description of the states is provided by the generalized state given in Sec. El 

We now test the behavior of the quasiholes in a generalized Laughlin state. This 
can be done as before, but now we must note that the gradient of the state will not only 
depend on the radial, but also on the angular position of the quasiholes. Furthermore, 
it will depend on the parameter /3 as defined in Eqs. (126]) and ( 12"7|) . which is used to 
improve the overlap. As shown in Fig. [12} for given parameters rj and Qo the same value 
for (3 optimizes simultaneously the ground state, the quasihole state, and the state with 
two quasiholes. We define 

fp(R,<j>) = (VcMld^cM) ■ (34) 

This function f@ is quite lengthy, so we expand it in R and give only the lowest term 
(0(R 2 )): 



8115904 + 2799526/3 2 - 7102^94958 (1 - (3 2 )f3 cos(20) 2 
* 17142924 + 4477401/3 2 R (35) 

From the expression we see that for a fixed and small value of R, fp oscillates around 
R 2 /2, such that the angular integration J Q n fp(R, 4>)d(p again will yield a Berry phase 
close to the encircled area, thus half of the Berry phase accumulated by a normal particle. 
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Figure 13. The effective charges q y (blue dotted) and q x (red dashed), and q e s = 
(<Z.t +<Zy)/2 (green solid) of quasiholes in the generalized Laughlin state as a function of 
the admixture j3 of higher angular momentum to the Laughlin state, for N = 4 (left) 
and N = 6 (right). 



Formally, we can capture this oscillating behavior by defining two effective charges 
q x and q y , depending on the direction in which the quasihole moves. With this and 
by generalizing to arbitrary loops the Berry phase in the limit of small radial positions 
R < 1, the Berry phase can be written as, 

7c = j> ^r(q x y, -q y x) ■ dr = (q x + q y )^, (36) 

where A is the encircled area and Stokes' theorem has been applied. An effective charge 
can be defined simply as, q c s = (q x + q y )/2. The values of q x ,q y , and g e fr are plotted 
as a function of /3 in Fig. [13] for different N. As can be seen, in all cases the effective 
charges are close to 1/2. For small values of /3, which represent realistic states of the 
system, the value of the charges gets closer to 1/2 as the number of atoms in the system 
is increased. In summary, the average charge q e g has only a minor dependence on (3, 
which decreases as N increases. Though not realized in our system, we note that in 
the limit /3 — > 1, both charges q x and g y again coincide due to the recovered cylindrical 
symmetry of the state \&ci. 

Finally we introduce two quasiholes into the generalized Laughlin state. Following 
a procedure similar to the one for the adiabatic case presented in the previous section, 
we extract the statistical phase angle of the quasiholes. The result as a function of /3 is 
shown in Fig. [TUfor N = 4 and closed paths of different radii. While the quasiholes in 
the bulk, R — 1, remain with an almost constant phase angle <yCstat ~ 0.51, the phase 
angle of quasiholes closer at the edge of the system have a stronger dependence on (3. 
Thus, we find that the presence of a certain degree of non-adiabaticity, (3 < 0.7 does 
not spoil the presence of anyonic quasihole excitations of the Laughlin state. 
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Figure 14. (Color online) Statistical phase angle v? s tat for two quasiholes in generalized 
Laughlin state as a function of the variational parameter /3. The radial position of the 
quasiholes is at Aj_ (blue solid line), 1.4Aj_ (red dashed line), and 1.8Aj_ (green dotted 
line) . 



7. Evolution with N of the energy gap 

The typical scenario of an experiment that has as a goal the realization of a specific 
strongly correlated state, is to first prepare an easily attainable initial state. Then, one 
may follow adiabatically a route in parameter space which ends up in the final desired 
state [161 [20]. This final state is expected to be robust, with a mean-life time larger 
than the time necessary to perform the measurements. A crucial ingredient necessary 
for the success of such approach is to have energy gaps as large as possible over all the 
ground states involved along the route. 

In this section we concentrate on the study of the energy gap over the Laughlin 
state. We consider first the adiabatic case and then study the effect of the perturbation 
on the energy gap. The gap is the energy difference between the ground state and 
the lowest excitation in the thermodynamic limit. We take the thermodynamic limit 
by increasing N keeping the chemical potential constant during the process. In the 
homogeneous 2D case this corresponds to fixing gN = constant, our calculations are 
performed fulfilling this relation. We will analyze the behavior of the gap for increasing 
N by means of our exact diagonalization calculations up to N = 7. 

For the Laughlin state in the adiabatic/symmetric shown in Fig. [2] for 

N = 4 and L — 12 there are two lowest excitations with L = 12 — 4 = 8 and L = 13 
depending on the value of rj. This is a general result for any N, the excitations have 
L = N(N — 1) — N and L = N(N — 1) + 1. The excitation with L = 13 is an excitation 
of the center of mass of the system. This is due to the incompressibility of the Laughlin 
state. Namely, as shown in Fig. [31 the state can increase its angular momentum without 
changing its interaction energy. We ignore this excitation, since we are interested in bulk 
excitations. This means that the linear left branch in the Laughlin region in Fig. [2] must 
be extended up to r) = 1 (where one has the largest gap). All the branches on the right 
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Figure 15. (Color online) (a) A/g, in units of Hu>±, of the Laughlin state as a function 
of N in the adiabatic/symmetric case, (b) Value of r) 2 computed for TV = 3,4, 5, 6, and 
7 in the adiabatic/symmetric case, compared to the prediction explained in the text, 
r/2 = l- 0.1(gN/N 2 ). 



that would lie below this line, are edge excitations of different polarity. In addition, the 
energy of the first excited state over the Laughlin denoted as A, with L = N(N — 1) 
(i.e., in the same L-subspace) coincides with this largest gap at rj — 1 where effectively 
there is no trap. 

As a consequence, to see the evolution of the gap over the Laughlin state it 
is sufficient to calculate the first two eigenvalues of the energy spectrum in the 
L = N(N — 1) subspace for each N. In Fig. [ToTa) we show the evolution of A/g with 
N. The tendency up to N = 7 is to asymptotically recover the value of 0.1 previously 
obtained by Regnault et al. [36] assuming a toroidal geometry and later reproduced 
by Roncaglia et al. |19j . Since we have taken gN = 6, then A must compensate this 
iV-dependence and consequently tend to zero as A ~ 1/N. In effect, our results imply 
that the bosonic Laughlin state is observable only for few number of particles. 

For possible practical implementations, it is also important to quantify the size of 
the parameter region where the Laughlin can be produced. Therefore, a good estimate 
of 7/2 can be obtained taking advantage of the abovementioned coincidence. r\i is the 
critical value where the energies of the Laughlin and L = N(N — 1) — N state cross, or 

(1 - r/ 2 )(L - N)twj ± + V 1 = (l- V 2)L niu ± + V (37) 

where L = N(N — 1) and V and V\ are E int (L ) and E int (L — N) (see Fig. [3]), 
respectively. Or, 

^ l -wk (v ° =a) - (38) 

In addition, V\ coincides with the energy difference A between the ground state and the 
first excitation in the Lq subspace which tends, as iV increases (see Fig. fl5]) (a). to 

A ~ O.l^fajj. = 0.1(6/N)hu± (39) 

and then 

qN 

V2 = l-A/(Nhw ± ) -1-0.1^2 . (40) 
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In Fig. IToT b) we compare the prediction of this formula and the computed r] 2 for different 
N. As can be seen the formula agrees very well with the numerically obtained values. 

Finally, Fig. [16] shows the change of the gap with decreasing Qq. There are some 
important differences between the symmetric and the perturbed cases. The perturbation 
mixes a large number of subspaces and now it is not possible to ignore the right branch. 
The initial (772) and final frequencies at the boundaries of the Laughlin region are shifted 
to smaller values. The largest gap, the one at the upper vertex of the triangle is nearly 
constant. At 772 the perturbation opens a gap where degeneracy occurs in the symmetric 
case. As a consequence, we conclude that on the one hand, the perturbation favors the 
observability, and on the other hand, the detection has to be restricted to a small number 
of particles. 



0.06 




Figure 16. (Color online) Evolution of the energy gap, in units of haj± with decreasing 
fio in the Laughlin region. 

8. Summary and conclusions 

We have studied the possibility of producing relevant strongly correlated quantum states 
as ground states of a system of ultracold two-level atoms subjected to an artificial 
gauge field. The focus is on the formation of Moore- Read (Pfaffian), Laughlin and 
Laughlin-quasiparticle states. Considering a few number of atoms we have shown 
by exact diagonalization methods, that large squared overlaps between Pfaffian-like, 
Laughlin-quasiparticle-like and Laughlin-like states and the ground state of the system 
are found even for fairly small values of the external laser intensity KIq/Er. Reducing 
the laser intensity, correspondingly reducing the Rabi frequency, induces deformation 
in the system, increasing the angular momentum of the states. The structure of these 
deformed states is such that the main properties of the original undeformed states is 
retained in a broad region of parameter space. An analytical representation of the 
ground states consisting on the original states supplemented by a term affected by 
Jastrow factors increasing in two units their angular momentum provides large overlap 
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with the numerical results, and allows to get analytical insights into the fractional 
behavior of the quasihole excitations. 

We have checked that quasihole states on the Laughlin and generalized Laughlin 
states can be produced by means of additional laser beams. We have studied the 
fractional charge and anyonic statistics of such quasiholes making use of the analytical 
representation of the states. Both, the effective charge and the statistical phase angle 
are close to the expected value of 1/2, even for the small number of particles considered 
here. Such fractional behavior is found whenever quasiholes are able to evolve in the 
bulk of the system, thus, the agreement with the expected behavior gets better as 
the size of the system is increased. The admixture of higher angular momentum in 
the generalized Laughlin state does not modify the fractional behavior of its quasihole 
excitation. However, as the creation of two quasiholes requires a higher artificial field 
strength t] than the one necessary to get into the (generalized) Laughlin regime, our 
analysis is relevant mostly in the weakly perturbed regime, whereas in the highly 
perturbed regime an exact numerical evolution of the adiabatic movement would be 
needed, falling beyond the scope of the present work. 

Concerning the observability of Laughlin-like states, we find that decreasing 
hQo/Eji shifts the spectrum to smaller values of rj, which are thus further away from 
the instability region 77 = 1, favoring the experimental conditions. In addition, keeping 
the chemical potential constant, i.e. Ng = constant, we find that the observability of 
the Laughlin state is reduced to small number of particles due to the fast decrease of 
the largest possible gap over the ground state as N increases. 

An interesting aspect for future investigation are the excitations in the Pfaffian-like 
regime, which might obey non-Abelian statistics. A signature for this behavior would 
be a degeneracy in the state with four quasiholes [U]. However, our present system is 
too small to test this physics. 
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Appendix A. Effective Hamiltonian 

Let us consider a system described by a single particle Hamiltonian given by H = Hq + V 
where H is solvable and V can be treated as a perturbation. In addition, let us assume 
that the eigenvalues of H Q are grouped in manifolds well-separated in anergy, i.e., the 
eigenenergies E ia have the following property, 

\E ia — Ej a \ <C \E ia — Ejp\ (A.l) 

where a and /3 denote different manifolds and the index % denotes different states inside 
a manifold. In our case we have two manifolds, the lowest energy one [a = 2) described 
by H22 and the most energetic one described by Hu {(3 = 1), both together play the 
role of H , 

H^( H - °) (A,) 



V H 22/ 
acting on the spinor shown below Eq.( [H]). 

The difference in energy, hQ , is much larger than the typical expected values of 
H 2 2 which are of the order of the recoil energy E R = h 2 k 2 /2M). Within the lowest 
manifold we make the LLL assumption. This scenario is valid in the case where the two 
internal states are uncoupled and the system evolves always in or ^2} • 

Instead, if there is coupling between the manifolds, the Hamiltonian is represented 
by H + XV where A is a dimensionless parameter and V, in general, has non-zero 
matrix elements inside each manifold as well as between them. As long as A is small, 
the structure of the well-separated manifolds and their degeneracy is preserved with 
slight modifications. Physically, the coupling between the two manifolds means that 
the motion of an atom in ^2) is slightly modified by a sudden and short time period in 
the other manifold. The high frequency dynamics driven by (E ia — E^/h (a 7^ (3) is 
averaged by the slow dynamics driven by (Ei a — Ej a )/h. In a way, the wave functions 
of the manifold-2 are "influenced" (or dressed) by the wave functions of manifold-1. In 
our case, the structure of V is, 



,#21 

with no diagonal elements. 

The main goal will be to obtain an effective Hermitian Hamiltonian H' that acts 
only on the unperturbed manifold-2, though having the same eigenvalues as those of H, 
and with zero matrix elements between the two unperturbed manifolds. In this way, we 
will be able to consider only the FD functios and ignore manifold-1. 

The Hermiticity and the coincidence on the eigenvalues with H are achieved if we 
consider a unitary transformation from H to H' as, 

H' = THT ] (A.4) 

where T = e tS and S = S\ Or 

H' = e iS He~ iS 

= H+ [iS, H} + ^ [iS, [iS, H]] + ... = H + \H[ + \ 2 H' 2 + ... (A.5) 
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where we have considered the expansion S = XS\ + A 2 ^ + . . . and the condition that 
S has zero matrix elements inside each manifold. Grouping the terms in orders of A, it 
is possible to solve Si step by step as functions of known quantities. One arrives to the 
expression: 

(ia\H'\ja) = E ia 5ij + ( ia | XV | jot ) + - (ia\XV \kj) {kj\XV \ja) 



2 



+ 



+ ... (A.6) 



_Ei a Efcry E!ja _ 

The first term represents the unperturbed levels in the manifold-2, the second one is 
the direct coupling between unperturbed levels in the manifold-2 and the third term is 
the contribution of the indirect coupling through the manifold- 1. In our case the last 
equation reduces to, 



( ia \H'\ja)^E l2 5 lj +^^(i2\V\kl)(kl\V\j2)(^-^-^ 
where we have approximated Ek\ by hQ and considered E i2 <C MIq or 



(A.7) 



\ia\H'\ja)= E l2 5 l3 - ^ < i2 \ V 2 \ j2 ) (A.8) 



o 



thus, 



H> = H 22 - (A.9) 
"o 

which is the results used in Eq. 012p in Section [2J The interaction term is considered 
part of the non-perturbed term together with the kinetic contribution, see Eq. ( fT3|) . 
The explicit form of the perturbation term H 2 \H\ 2 used is, 

__ „ / h 4 2x 2 h A k 2 x 2 h 4 k 4 x 2 h 4 

H 21 H i2 = I — — - - - rn - + r^rr^ a + 



AM 2 w 4 M 2 w e 16M 2 w 4 6AM 2 w 2 



ikxyh k y h 



4M 2 w 5 6AM 2 w 4 

ikxh 4 ik 3 xh 4 \ _ / xfi 4 ikyh 4 

dv + TTo 7 - TTTTo o d x 



4M 2 w 3 8M 2 w J y \M 2 w 4 8M 2 w 3 J x 

f k 2 h 4 k 2 x 2 h 4 \ 2 / h 4 x 2 h 4 \ 2 

+ V 4M 2 + 4M 2 w 2 ) y + V 4M 2 w 2 + 2M 2 w 4 ) x ^ ' 

one can show that this operator does not conserve L, as it connects L'-subspaces with 
V = L + A where A = 0, ±2, ±4. 
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